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We test different Big Bang Nucleosynthesis scenarios using the recent results on the Cosmic 
Microwave Background Anisotropies provided by the BOOMERanG and MAXIMA-1 experi- 
ments versus the observed abundances of *ffe, D and Li. The likelihood analysis, based on 
Bayesian approach, shows that, in the case of high deuterium abundance. Yd = (2.0±0.5)10~^, 
both standard and degenerate BBN are inconsistent with the CMBR measurement at more 
than 3(7. Assuming low deuterium abundance. Yd ~ (3.4±0.3)10~^, the standard BBN model 
is still inconsistent with present observations at 2a level, while the degenerate BBN results 
to be compatible. Unless systematics effects will be found in nuclide abundances and/or in 
CMBR data analysis this result may be a signal in favour of new physics like a large chemical 
potential of the relic neutrino-antineutrino background. 
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One of the main goals of modern cosmology is the 
knowledge of the energy density content of the uni- 
verse. The four parameters fls, ^cdm, and J^a, 
giving, respectively, the baryon, cold dark matter, 
neutrino and cosmological constant contributions to 
the total energy density, in unit of the critical density, 
and the Hubble constant, Hq = lOO/i Km s~^Mpc~'^ , 
enter as crucial parameters in several cosmological ob- 
servables. A well known example is provided by the 
baryon density parameter which plays an essential role 
in determining the abundances of light nuclides pro- 
duced in the early universe. An increasing interest 
has been also devoted to other issues, as structure for- 
mation and the anisotropy of the Cosmic Microwave 
Background Radiation (CMBR). As for our theoret- 
ical understanding of the Big Bang Nucleosynthesis 
(BBN), to test the theoretical models which describe 
these aspects of the hot Big Bang model, it is es- 
sential to have precise measurements of the several 
D,ih? . Furthermore, combining different cosmological 
observables, and comparing the way they are able to 
constrain the ^ih? , allows for a check of the consis- 
tency of our present understanding of the evolution of 
the universe. This can provide new hints on phenom- 
ena which took place at the macroscopic cosmologi- 
cal level, or rather related with the very microscopic 
structure of fundamental interactions. 

This letter represents a contribution in this direc- 
tion. In particular, we do perform a combined analysis 
of the dependence on the energy fractions Q^sh^ and 
Q,yh? = NyVl^h"^ for massless neutrinos {Ni, standing 
for the effective neutrino number and fi^ft.^ for the en- 
ergy contribution of a single v ~v specie) of CMBR 
anisotropies and BBN. This is aimed to test the stan- 
dard and degenerate BBN scenario, using the recent 
results of the BOOMERanG and MAXIMA-1 [§ 
CMBR experiments and the measurements of e, D 
and ''Li primordial abundances. 

The theoretical tools necessary to achieve this goal 



are nowadays rather robust and the accuracy in the 
predictions of the Big Bang model is remarkably im- 
proved. In fact, recently, a new generation of BBN 
codes have been developed which give the '^He 

mass fraction Yp with a theoretical error of few per 
mille. This effort is justified in view of the small sta- 
tistical error which is now quoted in the Yp measure- 
ments. On the other hand, the theoretical predictions 
on the CMBR anisotropies angular power spectrum, 
in the case of primordial passive and coherent per- 
turbations, expected in the most general inflationary 
model, have also recently reached a 1% level accuracy, 
with a refined treatment of H, He I, and He II recom- 
bination [0. 

Concerning the CMBR experimental data, after the 
COBE satellite first detection of CMBR anisotropies, 
at scales larger than 5°, and more than 20 indepen- 
dent detections at different frequencies and scales, see 
e.g. an important new insight is represented by the 
recent results obtained by the BOOMERanG Collab- 
oration [Q. For the first time, in fact, multifrequency 
maps of the microwave background anisotropies were 
realized over a significant part of the sky, with ~ 
10' resolution and high signal to noise ratio. The 
anisotropy power spectrum, C^, was measured in a 
wide range of angular scales from multipole ^ ^ 50 
up to ^ ~ 600, with error bars of the order of 10%, 
showing a peak at Ipeak — (197±6) with an ampli- 
tude DT2Q0 = {69±8) fiK . While the presence of such 
peak, compatible with infiationary scenario, was al- 
ready suggested by previous measurements the 
absence of secondary peaks after £ > 300 with a flat 
spectrum with an amplitude of ~ AOfj^K up to ^ ~ 625 
was a new and unexpected result. This result obtained 
then an impressive confirmation by the MAXIMA-1 
1^ experiment up to £ ~ 800. 

As far as BBN is concerned, in the last few years 
many results have been also obtained on light element 
primordial abundances. The ^He mass fraction Yp, 
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has been measured with a 0.1% precision in two in- 
dependent surveys, from regression to zero metalhcity 
in Blue Compact Galaxies, giving a low value Yp^'^ = 

0.234±0.003 0, and a high one Y^''^ = 0.244±0.002 
[n]] , which are compatible at 2a level only, may be due 
to large systematic errors. As in |^], in our analysis we 
adopt the more conservative value Yp = 0.238±0.005. 
A similar controversy holds in D measurements, where 
observations in different Quasars Absorption line Sys- 
tems (QAS) lead to the incompatible results Y^'-* = 

(3.4±0.3)10-5 H, and F^^^ = (2.0±0.5) 10"'' |l|. 
We will perform our analysis for both low and high D 
data. Finally, the most recent estimate for Li primor- 
dial abundance, from the Spite plateau observed in the 
halo of POP 11 stars, gives Yv^^ = (1.73±0.21) 10"^" 
PH ]. The light nuclide yields strongly depend on the 
baryon matter content of the universe, f^s/i^. High 
values for this parameter result in a larger '^He mass 
fraction and a lower deuterium number density. In 
particular, assuming a standard BBN scenario, i.e. 
vanishing neutrino chemical potentials, the likelihood 
analysis gives, at 95% C.L. B, 



low D neh^ = .017±.003 
high D ngh^ = .007; 



.007 



1.7 <N^< 3.3 
2.3 <N^< 4.4 



(1) 



As already pointed out by several authors ]15|-|17[|, 
these values for flsh^, though in the correct order 
of magnitude, are however somehow smaller than the 
baryon fraction which more easily fit the CMBR data. 
In fact, as mentioned, while the narrow first peak 
around I 200 is a confirmation of the infiationary 
paradigma, a fiat universe with adiabatic perturba- 
tions, the lack of observation of a secondary peak at 
smaller scales raises new intriguing questions about 
the values of the cosmological parameters. In partic- 
ular, this result may be a signal in favour of a larger 
Vlsh^ ~ 0.03, since increasing the baryon fraction en- 
hances the odd peaks only. In this respect the mea- 
surement of the third peak at larger multipole mo- 
ments is extremely important. Though it is fair to 
say that a further analysis of the BOOMERanG and 
MAXIMA- 1 data, as well as new data from future ex- 
periment, are needed to clarify this issue, it is however 
timely to study how our theoretical understanding of 
BBN can be reconciled with the larger values of flsh^ 
suggested by CMBR data. 

In figure 1 we start addressing the problem quan- 
titatively, by reporting our theoretical predictions on 
the several abundances, together with the correspond- 
ing experimental values, and the CMBR 68% C.L. 
bound on Qsh^ — 0.0301q'qq4, obtained by our joint 
analysis of both BOOMERanG and MAXlMA-1 data 
(see below for the details). The BBN predictions are 
obtained in the standard scenario for N^, = 3. Unless 
systematics effects will be found in nuclide abundances 
and/or in CMBR data analysis, we see that there is 
a rather sensible disagreement between the values of 
the baryon energy density required by standard BBN 
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FIG. 1. The theoretical abundances for D, "^He, '^He 
and '^Li nuclides, for A^^ = 3, versus the experimental 
results, represented by the horizontal bands, are reported. 
The long vertical band represents the 68% C.L. interval 
for ngh'^ obtained by the BOOMERanG and MAXIMA-1 
data analysis. 

and by CMBR data. In particular, from the 68% C.L. 
CMBR result for ilsh^, we can derive the following 
constraints on the Yi, 



0.2487 < YiHe< 0.2536 
9.23-10"^ < Yd < 2.3240"^ 
7.92-10-^ < Y3H^< 9.85-10-^ 
5.45-10-1° < Y7Li< 1.25T0-^ 



(2) 



which are systematically different from the available 
measured values. We stress once again that our the- 
oretical analysis is accurate, at worst, at the level of 
few percentj]. 

It has already been stressed ||l^,^ that a simple way 
to improve the agreement of observed nuclide abun- 
dances with fis/i^ > 0.02 is to assume non vanish- 
ing neutrino chemical potentials at the BBN epoch, a 



*It is important to note that this result is obtained fol- 
lowing a Bayesian approach in the analysis. Assuming the 
standard BBN, the best fit model for the CMBR spectrum 
has a normalized ~ 50, over ~ 40 degrees of freedom. 
Being the slightly greater than one, a CMB analysis 
using the BBN prior as in is still perfectly consistent. 



2 



scenario already extensively studied in the past ||T^ . 
The effect of neutrino chemical potentials fXa, with 
a the neutrino specie, is twofold. A non-vanishing 
Ca = Mi^q/^i'i contribute to N^, as 



implying a larger expansion rate of the universe with 
respect to the non-degenerate scenario, and a higher 
value for the neutron to proton density ratio at the 
freeze-out. Furthermore, a positive (negative) value 
for means a larger (smaller) number of Ve with re- 
spect to t'e, thus enhancing (lowering) n ^ p pro- 
cesses. Notice that extra relativistic degrees of free- 
dom, like light sterile neutrinos, would contribute to 
iVy as well, and in this respect BBN cannot distin- 
guish between their contribution to the total universe 
expansion rate and the one due to neutrino degener- 
acy. Therefore our estimates for N^, can only represent 
an upper bound for the total neutrino chemical poten- 
tials. However, it is also worth stressing that a large 
> 3, does indeed require a positive ^e, since to end 
up with the observed values for the 1^, a fine tuned 
balance between and iV^ is required [p"8P]. 

Increasing also weakly affects the CMBR 

anisotropy spectrum in two ways. The growth of per- 
turbations inside the horizon is in fact lowered, re- 
sulting in a decay of the gravitational potential and 
hence in an increase of the anisotropy near the first 
peak. Moreover, the size of horizon and sound hori- 
zon at the last scattering surface is changed, and this, 
with additional effects in the damping, varies the am- 
plitude and position of the other peaks, see e.g. p9[ . 

To test the degenerate BBN scenario we have per- 
formed a likelihood analysis of the data. First, to con- 
strain the values of the parameter set (^e, N,y, ^sh^) 
from the data on ^He, D and Li we define a total 
likelihood function^ 

jO-Nucl {N„,nBh'^ ,S,e) = Ld LiHe Ly Li , (4) 

as described in Ref. Q. For a fully degenerate BBN, 
since the effect of a positive can be compensated 
by larger N^, CNuci{N^,nBh'^,^e) may sensibly dif- 
fer from zero in a region with rather large values of 
N,y. We have chosen to constrain this parameter to 
be Ni, < 16. This upper limit has been considered 
after checking that it is well outside the 95% upper 
limit on from the BOOMERanG and MAXIMA-1 
data (again, see below). The other two parameters 
are chosen in the following ranges, — 1 < < 1 and 
0.004 < f^s/i^ < 0.110. 

Since CMBR spectrum is not sensible to alone, 
we have marginalized over ^g. All nuclide abun- 
dances have been evaluated using the new BBN code 
described in M, while the theoretical uncertainties 
al^'^{Nv, ^shr) are found by linear propagation of 
the errors affecting the various nuclear rates entering 
in the nucleosynthesis reaction network Q . 



The CMBR data analysis methods have been al- 
ready extensively described in various papers |p^ . 
The anisotropy power spectrum, was estimated in 12 
orthogonalized (independent) bins between £ = 50 
and £ = 650 by the BOOMERanG experiment and 
in 10 bins, from £ = 70 to ^ = 750, by the MAXIMA- 
1 experiment. The likelihood function of the CMBR 
signal Cb inside the bins is well approximated by 
an offset lognormal distribution, such that the quan- 
tity Db — ln{CB + xb) (where xb is the offset cor- 
rection) is a gaussian variable The likelihood 
for a given cosmological model is then defined by 
-21nL = (i?f - D%^)Mbb' {D%\ ~ D'j^,), where Mbb' 
is a matrix that defines the noise correlations and D^g 
is the offset lognormal theoretical band power. The 
theoretical were generated using fast and accu- 
rate Boltzmann solvers [|T|. Our database of models 
is sampled in physical variables as in [Ts] ] , but we only 
consider flat models, h — 0.65±0.2 and the effective 
number of neutrinos is allowed to vary up to N,y — 20. 
Following we can therefore constrain the param- 
eter of interest, U,Bh^ and N^, by finding the remain- 
ing "nuisance" parameters which maximize them. As 
mentioned, at 95 and 99% C.L. we found < 13 
and Ni, < 16, respectively. A similar bound, using 
BOOMERanG data only, has been obtained in 

In figure 2 we have summarized the second main re- 
sult of our analysis. In the Q,Bh^ — N^, plane we show 
the 95% C.L. likelihood regions for both the high and 
low D measurements, as well as the analogous con- 
tours for standard BBN, obtained running our code 
with = 0. In the same plot we show the 68 and 95 
% C.L. regions obtained by CMBR data. 

As a first comment, the standard BBN, = 0, and 
CMBR data analysis lead to quite different values for 
fisft,^. This can be clearly seen from the reported 95% 
results, but we have verified that the 99% C.L. con- 
tour for high D has no overlap with the region picked 
up by BOOMERanG and MAXIMA-1 data, and a 
very marginal one for low D. For the degenerate sce- 
nario, increasing N^, the allowed intervals for flBh^ 
shift towards larger values. However the high D val- 
ues require a baryon content of the universe energy 
density which is still too low, VlBh^ < 0.018, to be 
in agreement with CMBR results. A large overlap is 
instead obtained for the low D case, whose preferred 
^Bh"^ span the range 0.012 < VLbH^ < 0.036. As ex- 
pected, a larger helps in improving the agreement 
with the high CMBR fisft,^ value, but is important 
to stress that a large value for N^, is not preferred by 
the CMBR data alone, being, in this case, the best fit 
Ni, ~ 3. If we only consider the 95% overlap region 
we get the following conservative bounds: 

4 < < 13 , 0.024 < riBh'^ < 0.034 . (5) 

In this region varies in the range 0.07 < < 0.43. 
As we said, values iV^ > 3, as suggested from our anal- 
ysis, can be either due to weak interacting neutrino 
degeneracy, or rather to other unknown relativistic 
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FIG. 2. The 95% C.L. contours in the Qsh^-N^ plane 
compatible with degenerate BBN (large bands) and stan- 
dard BBN (small regions) are plotted for both high D (left) 
and low D (right). The same contours for 68% and 95% 
C.L. from BOOMERanG and MAXIMA-1 CMBR data 
are also reported. 



degrees of freedom. 

In conclusion, we have shown how a precision anal- 
ysis of BBN and CMBR data is able to tell us about 
possible new features of the cosmological model de- 
scribing the evolution of the universe. We have quan- 
titatively discussed how larger values for the bary- 
onic matter content Qsh^ may be reconciled with 
BBN predictions in a degenerate scenario. In this re- 
spect the observed low value of deuterium more easily 
fits with the constraints given by BOOMERanG and 
MAXIMA-1 data. The crucial aspect of our analy- 
sis is that this agreement is realized with an effec- 
tive neutrino degrees of freedom larger than 4 at 95%. 
This represents an indication in favour of a degener- 
ate neutrino background and/or new particle species 
contributing to relativistic matter in the universe. As 
final remark, we should note that our CMB analysis 
was restricted on a specific class of models with a lim- 
ited numbers of parameters. Including curvature, a 
gravity waves background, or removing our priors on 
h would not move the peak of our likelihood on Q,Bh^ , 
but would enlarge our C.L. more on the high value side 
. Including drastical deviations from the standard 
model like topological defects , decoherence in the 
primordial fluctuations [ p5[ or assuming a less general 
inflationary model would drastically change the 
conclusions of our work. Each of these effect leaves a 
characteristic imprint on CMB, so hopefully with new 
data available in the near future, as well as further 
analysis of the BOOMERanG and MAXIMA results, 
it will be possible to severely scrutinize this result. 

One of the author, A. Melchiorri would like to thank 
P. de Bernardis, P. Ferreira, J. Silk and N. Vittorio for 
valuable comments and discussions; G. Micle would 
like to thank CERN TH-Division for support. 



[1 
[2 

[3: 

[4" 

[5 

[6: 

[7 

[s: 

[9 

[lo; 
[11 

[12 
[13 

[14 
[15 



[16: 
[17: 



[is: 

[19 
[20 

[21 



[22: 



[23' 
[24' 
[25^ 
[26' 



P. de Bernardis et al. (Boomerang Coll.), Nature 404 
955 (2000). 

S. Hanany et al. (Maxima Coll.), Astrophys. J. , sub- 
mitted, (2000). 

S. Esposito, G. Mangano, G. Miele, and O. Pisanti, 
Nucl. Phys. B540 3 (1999); S. Esposito, G. Mangano, 
G. Miele, and O. Pisanti, Nucl. Phys. B568 421 
(2000). 

S. Esposito, G. Mangano, G. Miele, and O. Pisanti, 



istro-ph/0005571 



R.E. Lopez and M.S. Turner, Phys. Rev. D59 103502 

(1999) . 

E. Lisi, S. Sarkar, and F.L. Villante, Phys. Rev. D59 
123520 (1999). 

S. Seager, D.D. Sasselov, and D. Scott, Astrophys. 
J. 523 1 (1999). 

E. Pierpaoh, D. Scott, and M. White, Science 2171 

(2000) ; L. Knox and L. Page, Phys. Rev. Lett, in press, 
(2000). 

P. D. Mauskopf et al, Astrophys. J. , L59 536, (2000); 

A. D. Miller et al., Astrophys. J. , LI 524, (1999). 

K.A. Olive and G. Steigman, Astrophys. J. Suppl., 97 

49 (1995), and references therein. 

Y.I. Izotov and T.X. Thuan, Astrophys. J. , 500 188 

(1998). 

S. Buries and D. Tytler, Astrophys. J. , 499 689 
(1998). 

For the most recent analysis see J.K. Webb, R.F. 
Carswell, K.M. Lanzetta, R. Ferlet, M. Lemoine, A. 
Vidal-Madjar, and D.V. Bowen, Nature, 388 250 
(1997). 

P. Bonifacio and P. Molaro, Mon. Not. Roy. Astron. 
Soc, 285 847 (1997). 
A.E. Lang e et al. (Boomerang 



Coll.' 



astro- 



astro- 



ph/0005004 ; M. Tegmark and M. Zaldarriaga^ 
ph/0004393 , A. H. Jaffe et al., |a.stro-ph/0007333|; W 
Kinney, A. Melchiorri, A. Riotto, [a,stro-ph/000737£| . 



M. Le Dour, M. D ouspis, J.G. Barlett, A. Blanchard, 
astro-ph/0004283l M.White, D. Scott, and E. Pier- 



paoh, |astro-ph/0004385 



J. Lesgourgues and S. Pastor, Phys. Rev., D60 103521 
(1999); W.H. Kinney and A. Riotto, Phys. Re v. Lett 
83 3366 (19 99); J. Lesgourgues and M. Peloso, 



astro- 



ph/0004412; M. Orito, T. Kajino, G.J.Mathews, and 



R.N.Boyd, |astro-ph/0005446 . 

H. Kang and G. Steigman, Nucl. Phys. , B372 494 
(1992). 

W. Hu et al., Phys. Rev., D52 5498 (1995). 

J.R. Bond, A. H. Jaffe, and L. Knox, Astrophys. J. , 

533 19 (2000). 

U. Seljak and M. Zaldarriaga, Astrophys. J. , 469 437 
(1996); A. Lewis, A. Challinor, and A. Lasenby, ^s- 
trophys. J. , in press (2000). 

L. Knox and S. Dodelson, Phys. Rev. Lett., 84 3523 
(2000); A. Melchiorri et al., Astrophys. J., 536 L63 
(2000). 



S. Hannestad, |astro-ph/0005018 
F. R. Bouchet et al, |astro-ph/0005022 



A. Lewin, and A. A lbrecht, |astro-ph/9908061 
T. Kanazawa et al, ^stro-ph/000644 



4 



